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I ' Abstract 

^ ; 

^. We have calculated the superconducting critical temperature Tc of hole-doped 

c/2 ! delafossite CUAIO2 based on the first-principles calculations. According our 

calculation, 0.2 ~ 0.3 hole-doped CUAIO2 can become a phonon-mediated 
high-Tc superconductor with Tc ~ 50 K. In the hole-doped CUAIO2, the AiLi 
phonon mode that stretches 0-Cu-O dumbbell has a strong interaction with 
r^ • electrons of the flat band in Cu 3d322_r2 and the O 2p2 anti-bonding vr-band. 



1. Introduction 



Cj I Keywords: A. Semiconductors; C. Delafossite structure; D. 

Electron-phonon interactions; E. Density functional theory 

m 
> 

^ ' Kawazoe et al. have discovered the delafossite structure of CUAIO2 is the 

t^ ■ transparent p-type conductor without any intentional doping, jlj] Transparent 

O \ p-type conductors such as CUAIO2 are rare and absolutely necessary for the 

p-n junction of the transparent conductors and high-efficient photovoltaic 

solar cells. Many applications of CUAIO2 for flat panel displays, photovoltaic 
Lj solar-cells, touch panels, and high efficiency thermoelectric-power materials 

'^ ■ with about 1% hole-doping |2|, |3[ are expected. 

•^ ' Recently, Katayama-Yoshida et al. have suggested a new application of 

CUAIO2 for transparent superconductivity and high-efficient thermoelectric; 

power material with a large Seebeck coefficient caused by the fiat band. 
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They have simulated the p-type doped CUAIO2 by shifting the Fermi level 
rigidly with FLAPW method, and proposed that the nesting Fermi surface 
may cause a strong electron-phonon interaction and a transparent supercon- 
ductivity for visible light due to the large band gap (~ 3.0eV). But, the 
calculation of superconducting critical temperature Tc is not carried out. In 
this study, we calculated the electron-phonon interaction and the T^ of p- type 
doped CUAIO2 based on the first principles calculation with the density func- 
tional perturbation theory. J5| We found that the Tc goes up to about 50 K 
due to the strong electron-phonon interaction and high phonon frequency 
caused by the two dimensional fiat band in the top of the valence band. 

2. Calculation Methods 

The calculations are performed within the density functional theory|6|, |7| 
with a plane-wave pseudopotential method, as implemented in the Quantum- 
ESPRESSO code.jsj We employed the Perdew-Burke-Ernzerhof generalized 
gradient approximation (GG A) exchange-correlation functional [9| and ultra- 
soft pseudopotentials. [10l| For the pseudopotentials, Cu 3d electrons were 
also included in the valence electrons. In reciprocal lattice space integral 
calculation, we used 8x8x8 (electron and phonon) and 32 x 32 x 32 



(average at Fermi level) k-point grids in the Monkhorst-Pack grid.[ll|| The 
energy cut-off for wave function is 40 Ry and that for charge density is 320 
Ry. These k-point meshes are fine enough to achieve convergence within 0.1 
mRy/atom in the total energy. The 32 x 32 x 32 mesh for average at Fermi 
level is enough to achieve convergence in the electron-phonon interaction and 
the superconducting critical temperature. The differences between results of 
32 X 32 X 32 k-points mesh and those of 64 x 64 x 64 one are less than 1%. 
The delafossite structure belongs to the space group R3m (No. 166) and 
is represented by cell parameters a and c, and internal parameter z (See 
Fig. [T]). These cell parameters and internal parameter were optimized by 
the constant-pressure variable-cell relaxation using the Parrinello-Rahman 
method [12| without any symmetry requirements. The results of relaxation 
(a = 2.86lA,c/a = 5.969 and z = 0.1101) agree very well with the experi- 



mental data (a = 2.858A,c/a = 5.934 and z = 0.1099 [13|, Ijj). 

In this study, some properties of hole-doped CUAIO2 are approximated 
because it is difficult for first-principles calculation to deal with the doped 
system exactly. Let's take the electron-phonon interaction A for example. A 
is defined as follows: 
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Figure 1: The crystal structure of delafossite CuA102- 



(1) For the non-doped CUAIO2, we calculated the dynamical matrix, the 
phonon frequency u^^ and the electron-phonon matrix M^*^^ . (2) For the 
doped CUAIO2, we calculated the Fermi level ep and the density of states at 
the Fermi level N{e-F) with the number of valence electrons reduced using 
the eigenvalues ^k of the non-doped system. (3) By using the results of 
(1) and (2), we calculated the electron-phonon interaction A and the other 
superconducting properties. This approximation is based on the rigid band 
model and the idea that the doping does not greatly change the phonon 
band structures. In this study, we show the results of 0.1 ~ 1.0 hole-doped 
CuAlOs. 

3. Calculation Results and Discussion 

Before the discussion of the superconducting critical temperature, let us 
see the electronic structure. Fig. [2] and Fig. [3] show the electronic band 
structure and density of states (DOS) of non-doped CUAIO2. The top of 
the valence band of CUAIO2 is fiat due to the two dimensionality in 0-Cu-O 
dumbbell array, and has a small peak in the DOS of the valence band. This 
peak is mainly constructed by the two-dimensional 7r-band of Cu 3^322 _r2-0 
2pz anti-bonding state. 

Fig. m shows the DOS at the Fermi level calculated with the number of 
valence electron reduced. According to this figure and Fig. |3l the number of 
holes iVh = 0.3 corresponds to that Fermi level which is located at the top of 
the peak in the DOS, and N-^ = 0.9 corresponds to that Fermi level which is 
located at the bottom of the DOS. 

We calculated the superconducting critical temperature by using the 



Allen-Dynes modified McMillan's formula. |15l.ll6l| According to this formula, 
Tc is given by three parameters: the electron-phonon interaction A, the loga- 
rithmic averaged phonon frequency cuiog, and the screened Coulomb interac- 
tion fi*, in the following form. 

^^ " 1.2 """"^y A-/x*(l + 0.62A)y'- ^^^ 

wiog = exp I - / duj logo; I, (3) 



Here, a^F{u)) is the Eliashberg function. A and uj\og are obtained by the first- 
principle calculations using the density functional perturbation theory. As 
for /i*, we assume the value n* = 0.1. This value holds for weakly correlated 
materials due to the electronic structure of lightly hole-doped Cu+ (d^°). 

The calculated result of Tc and A as a function of A^h is shown in Fig. |5l 
In our calculation, the lightly doped CuAlOa {N^ = 0.2 ~ 0.3) has T^ ~ 50 K. 
This Tc is the highest among phonon-mediated superconductors. In addition, 
the Tc can be increased by other purely attractive electron-electron interac- 
tion mechanisms: for example, charge-excitation-induced ITJ or exchange- 



correlation-induced [l8| negative effective U system in the Cu+ (d^°) elec- 
tronic structure with lightly hole-doping. The heavily doped CUAIO2 (A^^ = 
0.6 ~ 1.0) has Tc ^ lOK by reducing the electron-phonon interaction. 



4 
2 


-2 
-4 
-6 
-8 
-10 






Figure 2: Band structure of CuAlOa 
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Figure 3: Total density of states(DOS) and projected DOS. 
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Figure 4: The number of holes vs density of states at the Fermi level. 
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Figure 5: Superconducting critical temperature and electron-phonon interaction A. 



0.3 0.931 789 
1.0 0.405 778 

Table 1 : Electron-phonon interaction A and logarithmic averaged phonon frequencies wiog 
Tc has max. and min. at iVh = 0.3, 1.0. 



Let us examine the origin of the high T^ of hghtly hole-doped CUAIO2. In 
this study, the critical temperature is determined by A and wiog as mentioned 
above. Table |3]shows A and wiog. The electron-phonon interaction A at A^h = 
0.3 is about 130% larger than that at A^h = 1-0, while uiog at A^h = 0.3 is about 
1% larger than that at A^h = 1-0. In addition, A affects Tc exponentially, 
while uiog affects Tc hnearly. Therefore, the high Tc is attributed to the 
strong electron-phonon interaction. 

In Fig. ini we show the phonon dispersion which show the strong two di- 
mensionality with flat phonon dispersion. In order to find that which phonon 
mode has a large contribution to the high Tc, we introduce a partial electron- 
phonon interaction X^^'- the interaction of a phonon whose frequency is Wj^q. 
Then, A = ^^^ Aj^q. In Fig. [6], the A,^q is shown by the radius of a circle on 
each phonon dispersion. Since most of X^^ are very small, their circles are no 
longer invisible in the figure El This figure indicates that the highest mode 
on the Z — r line has a large electron-phonon interaction. In the case of 
A^h = 0.3, the sum of Aj^q of the highest frequency mode is 0.407. This value 
is about 44% of total electron-phonon interaction A = 0.931. The effective 
mode is the AiLi phonon mode. In this mode, the O atoms oscillate in the 
anti-phase within an 0-Cu-O dumbbell. 

As mentioned above, the CUAIO2 has the Cu 3(i3^2_^2 and the O 2pz 
electrons at the top of the valence band. When A^h = 0.2 ~ 0.3, the electrons 
which make the 0-Cu-O anti-bonding band are located at the Fermi level. 
They have a strong interaction with the AiLi phonon mode because their 
bonding direction is parallel to the oscillation direction of the AiLi phonon 
mode. Though the strong electron-phonon interaction, the 0-Cu-O bonding 



of the delafossite structure is stable even under high pressure. [19[] There is 
a strong possibility that the doped CUAIO2 is stable and a superconductor. 
When CUAIO2 is heavily hole-doped, the total electron-phonon interaction 
decreases because the number of electrons which have a strong interaction 
decreases. 





(b) 



Figure 6: Phonon dispersions and electron-phonon interactions of hole-doped CuA102. 
The radius of circle represents the strength of partial electron-phonon interaction Aj/q. 
Note that many A^q are very small and their circles are no longer invisible, (a) The 
number of holes TVh = 1.0. (b) iVh = 0.3. 



The top of valence band is constructed by the Cu 3^322 .^2 and the O 2pz 



anti-bonding vr-band. [19|, l20|, |2l| The hole-doping makes the 0-Cu-O coupling 



more strong. Therefore, when the density of holes increases from 0.2 to 0.3 
A and wiog does not change much (A = 0.901 — )■ 0.931, wiog = 808 — > 789K). 

D. Huang and Y. Pan have investigated the intrinsic defects in CUAIO2. 22 
According to their study, vacancies at the Cu sites and substitutional Cu at 
the Al site are most likely responsible for the p-type conductivity, and the 
transition levels of these defects are deep. The rigid-band doping is possibly 
not realized in CUAIO2. Many theoretical proposals for superconductivity- 
in-semiconductors (for example, LiBC 23|) are based on rigid band or other 
oversimplified models and not successful in experiment because doping con- 
centrations cannot be realized or large doping levels cause material distor- 
tions not accounted for by these models. Therefore, the calculated T^ — '^QK 
may not be realized in experiment. However, we believe our calculation re- 
sults suggest the superconductivity potential of hole-doped CUAIO2. 



4. Conclusions 



In summary, we calculated the superconducting critical temperature of 
the hole-doped delafossite CUAIO2 by shifting the Fermi level rigidly based 
on the first principles calculation. The lightly hole-doped CUAIO2 has Cu 
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■iz^ 



r2 and O 2pz anti-bonding vr-band as the top of the valence band. 



The electrons of this band have a strong electron-phonon interaction with 
the AiLi phonon mode because the direction of 0-Cu-O dumbbell is par- 
allel to the oscillation direction of the AiLi phonon mode. These findings 
suggest that hole-doped CUAIO2 may be a superconductor. We hope that 
our computational materials design of superconductivity will be verified by 
experiments very soon. We can easily extend the present computational ma- 
terials design to other delafossite structures of lightly hole-doped AgA102 
and AUAIO2 which may have higher Tc due to the strong electron-phonon 
interaction combined with the charge excitation- induced [l7J] and exchange- 



correlation- induced [l8| negative effective U; such as 2Ag^+ (d^) — )■ Ag+ (d^°) 
+ Ag3+ (d^) and 2Au2+ (d^) ^ Au+ (d^") + Au3+ (d^)^ upon the hole- 
doping. 

Acknowledgement 

The authors acknowledge the financial support from the Global Center 
of Excellence (COE) program "Core Research and Engineering of Advanced 
Materials - Interdisciplinary Education Center for Materials Science", the 
Ministry of Education, Culture, Sports, Science and Technology, Japan, and 
a Grant-in-Aid for Scientific Research on Innovative Areas "Materials De- 
sign through Computics: Correlation and Non-Equilibrium Dynamics". We 
also thank to the financial support from the Advanced Low Carbon Technol- 
ogy Research and Development Program, the Japan Science and Technology 
Agency for the financial support. 

References 

[1] H. Kawazoe, M. Yasukawa, H. Hyodo, M. Kurita, H. Yanagi, H. Hosono, 
Nature 389 (1997) 939. 

[2] H. Funashima, A. Yanase, H. Harima, H. Katayama-Yoshida, Proceed- 
ings of 23rd International Conference on Thermoelectrics (2004) 237- 
238. 

[3] H. Katayama-Yoshida, H. Funashima, I. Hamada, H. Harima, 
A. Yanase, Fabrication method and high-efficient thermo-electric power 
materials by two-dimensional natural super- lattices of cualo2, Japan 
Patent 652696 (2009). 



[4] H. Katayama-Yoshida, T. Koyanagi, H. Funashima, H. Harima, 
A. Yanase, Solid State Communication 126 (2003) 135. 

[5] S. Baroni, S. de Gironcoli, A. D. Corso, P. Giannozzi, Rev. Mod. Phys. 
731 (2001) 515. 

[6] P. Hohenberg, W. Kohn, Phys. Rev. 136 (1964) B864. 

[7] W. Kohn, L. J. Sham, Phys. Rev. 140 (1965) A1133. 

[8] P. Giannozzi, S. Baroni, N. Bonini, M. Calandra, R. Car, C. Cavazzoni, 
D. Ceresoh, G. L. Chiarotti, M. Cococcioni, I. Dabo, A. D. Corso, S. Fab- 
ris, G. Fratesi, S. de Gironcoh, R. Gebauer, U. Gerstmann, C. Gougous- 
sis, A. Kokalj, M. Lazzeri, L. Martin-Samos, N. Marzari, F. Mauri, 
R. Mazzarello, S. Paohni, A. Pasquarello, L. Paulatto, C. Sbraccia, 
S. Scandolo, G. Sclauzero, A. P. Seitsonen, A. Smogunov, P. Umari, 
R. M. Wentzcovitch, J. Phys.:Condens. Matter 21 (2009) 395502. 



[9 
[10 

[11 
[12 
[13 

[14 
[15 
[16 
[17 

[18 
[19 



J. P. Perdew, K. Burke, M. Ernzerhof, Phys. Rev. Lett. 77 (1996) 3865. 

D. Vanderbilt, Phys. Rev. B 41 (1990) 7892. 

H. J. Monkhorst, J. D. Pack, Phys. Rev. B 13 (1976) 5188. 

M. Parrinello, A. Rahman, Phys. Rev. Lett. 45 (1980) 1196. 

X. Nie, S.-H. Wei, S. B. Zhang, Phys. Rev. Lett. 88 (2002) 066405. 

A. Buljan, P. Alemany, E. Ruiz, J. Phys. Chem. B 103 (1999) 8060. 

W. L. McMillan, Phys. Rev. 167 (1968) 331. 

P. B. Allen, R. C. Dynes, Phys. Rev. B 12 (1975) 905. 

H. Katayama-Yoshida, K. Kusakabe, H. Kizaki, A. Nakanishi, Appl. 
Phys. Express 1 (2008) 081703. 

H. Katayama-Yoshida, A. Zunger, Phys. Rev. Lett. 55 (1985) 1618. 

A. Nakanishi, H. Katayama-Yoshida, J. Phys. Soc. Jpn. 80 (2011) 
024706. 



10 



[20] A. Nakanishi, H. Katayama-Yoshida, J. Phys. Soc. Jpn. 80 (2011) 
053706. 

[21] I. Hamada, H. Katayama-Yoshida, Physica B 376 (2006) 808. 

[22] D. Huang, Y. Pan, Can. J. Phys. 88 (2010) 927. 

[23] H. Rosner, A. Kitaigorodsky, W. E. Pickett, Phys. Rev. Lett 88 (2002) 
127001. 

[24] J.-Y. Son, T. Mizokawa, J. W. Quilty, K. Takubo, Phys. Rev. B 72 
(2005) 235105. 



11 



